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Abstract 



, We investigate a proposal for the construction of models with chiral fermions on the 

I \ lattice using staggered fermions. In this approach the gauge invariance is broken by the 
Q \ coupling of the staggered fermions to the gauge fields. Motivated by previous results in the 
■ non-gauge invariant massive Yang-Mills theory and certain gauge-fermion models we aim at 
a dynamical restoration of the gauge invariance in the full quantum model. If the gauge 
^ , symmetry breaking is not too severe, this procedure could lead in the continuum limit to 
^ ! the desired gauge invariant chiral gauge theory. This scenario is very attractive since it does 
not rely on gauge fixing. We investigate a simple realization of this approach in a U(l) 
axial-vector model with dynamical fermions in four dimensions. 
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1 Introduction 



An important unsolved problem in lattice field theory is the non-perturbative formulation of 
a chiral gauge theory on the lattice (for reviews see refs. fl], 3). The naive lattice transcription 
of a quantum field theory involving fermions leads to unwanted fermion species, the so-called 
fermion doublers. Half of these couple with opposite chiral charge to the gauge fields, spoiling 
the chiral nature of the couplings, and the resulting theory in the scaling region is vector- 
like. To deal with these doublers, one has three options: I) Remove the fermion doublers by 
rendering them heavy, while keeping only one fermion in the physical spectrum [3-11], II) 
decouple them by turning off their interactions with the other particles |jl2| or III) use them 
as physical degrees of freedom |13|, |l^ . 

When applied to chiral models, one has to reconcile the chosen regularization with gauge 
invariance. Many proposals are formulated such that local gauge invariance on the lattice 
is preserved [3,4,7-9,11]. For example, using a Wilson term to decouple the doublers, one 
introduces extra scalar fields to make this mass term gauge invariant. Alternatively one 
can sacrifice gauge invariance of the lattice model and transcribe the gauge fixed continuum 
action to the lattice [0, 0. This avoids introducing extra scalar fields, but requires Fadeev- 
Popov and gauge fixing terms. It can be done with Wilson fermions as in refs. [^, but 
also with staggered fermions. In either case one breaks the BRST symmetry, and one has 
to add counterterms to the action and tune their coefficients such that BRST invariance is 
restored in the continuum limit. This method is cumbersome from a technical point of view 
and in addition one has to worry about non-perturbative gauge fixing. 

The method which we will focus on in this paper falls in class III. It relies on stag- 
gered fermions, which method uses the species doublers as Dirac-flavor components. The 
one-component staggered fermion fields do not carry explicit Dirac and flavor labels, these 
components are 'spread out' over the lattice. In the classical continuum limit this method 
leads to four flavors of Dirac fermions. It is possible to couple these Dirac flavors to the 
gauge fields such that the chiral target model is recovered in the classical continuum limit. 
We studied the method in ref. for the case of a two-dimensional model with axial-vector 
couplings and concluded that it performed well for smooth external gauge fields. The im- 
portant issue, however, is whether the same can be achieved also when the full quantum 
fluctuations are taken into account. The same method of coupling the staggered Dirac-flavor 
components has recently been successfully used for an investigation of a strongly coupled 
fermion- Higgs model |TB[ . 

In chiral gauge models we are confronted with the problem that the couplings of the 
staggered fermions to the gauge fields spoil the local gauge invariance. One possibility to 
overcome this difficulty is described in ref. [Q, using the full machinery of gauge fixing. 
Alternatively one can attempt to avoid this and hope for a dynamical restoration of gauge 
invariance in the model without gauge fixing. By this we mean the following: carrying out 
the integration over all gauge field configurations in the path integral, the gauge degrees 
of freedom appear in the action as dynamical scalar fields with frozen radial mode; when 
these scalar fields decouple, gauge invariance is restored. As we shall see later there exist 
examples of non-chiral models where these scalar fields indeed decouple for suitable choices 
of the bare parameters. It is the subject of this paper to investigate this scenario for the 
staggered fermion approach to chiral gauge theories. 

The outline of the paper is as follows: In sect. 2 we recall the scenario of dynamical 
gauge symmetry restoration and review the massive Yang-Mills model and a gauge-fermion 
model, as examples of models in which this scenario was found to work. In sect. 3 we define 
two staggered fermion models with axial-vector couplings to a U(l) gauge field: the ISF 



(invariant staggered fermion) model, which has a local gauge invariance and the NISF (non- 
invariant staggered fermion) model, which couples the staggered Dirac-flavor components in 
a way that allows for the construction of chiral models, but breaks local gauge invariance. 
The phase diagram of the ISF model will be presented in sect. 4. We shall also break the 
gauge symmetry of the ISF model by hand and argue that in this case the gauge invariance 
can again be restored dynamically. The phase diagram of the NISF model is presented in 
sect. 5 and from it we infer that the desired symmetry restoration does most probably not 
take place. We try to improve on this by modifying the model in sect. 6. Also here we 
find no evidence for the dynamical restoration of the gauge symmetry. Sect. 7 contains a 
summary of our results and gives an outlook to possible future investigations. 



There exist examples of gauge non-invariant lattice models where a dynamical restoration 
of gauge symmetry takes place. More precisely it has been found that one can add terms 
to a gauge invariant action that break the gauge symmetry, but provided the bare coupling 
constants of these terms arc not too large, the low energy model remains the same. In our 
case the gauge symmetry is broken by the lattice regularization with staggered fermions, and 
there is no easy separation between a gauge invariant action and additional gauge symmetry 
breaking terms. The question arises if this symmetry breaking is sufficiently small that the 
chiral symmetry gets restored dynamically. 

To explain the basic idea let us start first from a generic non-gauge invariant lattice 
action Seff{U'). If also fermions are involved Seff{U') represents the effective action after 
integrating out the fermionic degrees of freedom in the path integral. The link field Uj^x — 
exp{—iaA^x) G G, with a the lattice distance, Aij,^ the vector potential and G the compact 
gauge group. In the following we shall mostly use lattice units, a — 1. Assuming that the 
partition function is defined by the usual integration over the link field, it can be written as 



where we have used in the second equation the gauge transformation U'^^ = Q;j:Ufjj:Qlj^_p as 
a transformation of variables as well as the invariance of the Haar measure, DU' = DU. In 
the third equation we have added the trivial integration 1 = / Dfl, and wrote Vx = ^l- In 
this way the gauge degrees of freedom Q have been exhibited as an additional field in the 
path integral, which could be interpreted as a dynamical Higgs fields V E G with a frozen 
radial mode, V^V^ = 1. The new action S'e//(V^^t/^j;\4+/i) is invariant under the local gauge 
transformations f/^^ QxU^x^l^^,, 14 ^xV^. 

The above manipulations turn any gauge non-invariant model into a related gauge invari- 
ant model, at the price of introducing an extra scalar field & G. The important question 
is if the resulting model still describes the physics of the underlying gauge invariant model, 
i.e. the models without adding the symmetry breaking terms to the action or, in our case, 
the model of the chirally invariant target model in the continuum. For instance, this requires 
that the V field decouples in the scaling region. 

Before we turn to the question of dynamical gauge symmetry restoration in the staggered 
fermion model, let us first present two specific examples, namely the massive Yang-Mills 
theory and certain gauge-fermion models, which show that gauge non-invariant terms with 



2 Restoration of gauge symmetry 




(2.1) 



coefficients which are not too large, do not spoil gauge invariance after the integration over 
quantum fluctuations in the path integral has been carried out. 

2.1 Massive Yang-Mills model 

The massive SU{2) Yang-Mills model is defined on the lattice by the euclidean action, 

Su = ^E'^rUx,. + Kj2^r{u,x + Ul,} , (2.2) 

where g is the gauge coupling and k, the bare mass parameter in lattice units (a = 1). The 
Uxfiu denotes the usual plaquette variable on the lattice. The mass term in ( |2.2|) breaks the 
local SU(2) gauge invariance. After going through the steps of eq. (|2.1| ) we find the action 
of the gauge invariant extension of the model, 

Suy = 4 E Tr f/x;.. + E Tr {v^U,,V,,^f, + V^^^UlV,} , (2.3) 

9 XflU Xfl 

which is the action of a gauged SU(2) xSU(2) non-linear sigma model. Note that the per- 
turbative non-renormalizability of the model ( |2.2| ) causes no problems when the model is 
treated non-perturbatively. 

The non-linear sigma model at g = (i.e. for f/^^. = 1) has a phase transition at 
K = Kc(0) ~ 0.3. For g > this phase transition extends into the n-g plane and separates a 
Higgs {k > Kc{g)) from a confinement phase (k < Kc{g)). We are interested here only in the 
scaling region at small bare gauge coupling g, where we find continuum behavior. 

There are three scaling regions of interest: (A) When approaching the point Kc(0) from 
within the Higgs phase, we encounter the usual Higgs phenomenon. The spectrum contains 
the three massive gauge bosons and the Higgs particle. (B) When approaching the point 
Kc(0) from within the confinement phase we have a theory with confined scalar particles, 
much like QCD. (C) If we let g approach zero inside the confinement phase away from the 
Higgs- confinement phase transition ndg), the scalar particles acquire masses of the order 
of the cut-off and we recover the pure SU(2) Yang-Mills system with only glueballs in the 
physical particle spectrum. 

In case C there is a whole region in parameter space (0 < k < Kc{g)), where the physics 
in the scaling region is the same as at k = 0, at which point the bare action is that of the 
gauge invariant Yang-Mills model. However, if we make the symmetry breaking too large 
(i.e. K close to Kc), the model changes and goes into a different phase in which the physics 
is no longer described by the pure Yang- Mills part of (|2.2|). It is also interesting to note 
in passing that in case (A) the gauge degrees of freedom V produce a genuine Higgs field, 
which acquires a non-frozen radial mode in the low energy action. The representation of the 
gauge group carried by this Higgs field is determined by the initial symmetry breaking. 

2.2 Gauge-fermion models 

As a second example we discuss restoration of gauge invariance in models with fermions in 
which the symmetry breaking resides in the fermionic part of the action. The model which 
we shall consider here, and in more detail in the next sections, has axial- vector couplings to 
U(l) gauge fields, 

S = - d'xY. {^'Md, + iqn,A^Wf + . (2.4) 

/=i 



There are Np fermion flavors with charges qj = ±1, / = 1, . . . ,Nf, such that the model 
is anomaly free. For y = the action (|2.4|) is invariant under the local gauge transforma- 
tion: A^{x) ^ A^{x) + d^tuix), ^(x) ^ in{x)PL + n*{x)PR)^ix), ^(x) ^(x)(fi(x)PL + 
Q*{x)Pji), with Q{x) = expiuj{x) and = (1 — 75)/2, Pr = (1 + 75)/2. A non-zero mass 
parameter |/ allows us to discuss the effect of gauge symmetry breaking by a fermion mass. 
We have put a prime on ip and ip to indicate the gauge in which the symmetry breaking 
has the simple mass term form J2fy4'f'^'f- Our target model is the gauge invariant model 
at y = 0. We shall show that this can be achieved even for nonzero y, provided it is not too 
large, 

The lattice transcription of the action ( p.4|) with a single naive fermion field reads 



Xfl 

-yj:^'.€- (2.5) 



Because of the species doubling phenomenon this action reduces in the classical continuum 
limit to ( |2.4| ) with Np = 16, where eight of the fermion species couple to the gauge fields 
with Qf = +1 and the remaining eight with qf = —1. Note that with these charges the 
model is anomaly free and equivalent to QED with 16 flavors. The model we shall actually 
study numerically is a staggered fermion reduction of ( p.5|) in which the number of flavors 
is reduced by a factor of two (four q/s equal to -|-1 and four equal to —1). This model we 
shall call in the following the ISF model. However, in order not to overload the reader here 
with the details of the staggered fermion formalism, we shall continue for the moment with 
( pl5|) as its properties are qualitatively the same as the ISF model. 

For y = the action (|2.5| ) is invariant under the local gauge transformations. 



U,. ^ Q^U^Ml^f, , ^ {Q^Pl + Q:Pr)^: , ^',-^i,'A^.PL + KPR) ■ (2.6) 

For non-zero y the mass term in the action breaks gauge invariance. We shall argue however, 
that as in the case of the massive Yang-Mills model, the physics in the scaling region remains 
the same as with ?/ = 0, provided y is not too large. Our arguments below are based on the 



knowledge of Wilson- Yukawa models accumulated in recent years [jT7|. 



Starting from the gauge non-invariant action ( |2.5| ) and after going through similar steps 
as in ( |2.1| ) we find. 



-i^UMv:+f.u;,v.PL + v:u,,v.+^PnW,} - • (2.7) 



The fields -ip' and tp' are now seen to be screened from the gauge fields by the V fields and 
therefore are neutral with respect to the U(l) gauge transformations 

u^^ n^u^^ni^^ , 14 ^xV^ , ^ , ^1 ^ i^L- (2.8) 

To see how the symmetry gets restored it is useful to make a transformation of variables to 
new fermionic fields 



V'^ = (vTPl + i4Pij)^. , ^1 = ^.(v;*Pl + v;Pr), 



(2.9) 



which leaves the integration measure invariant, DipDtjj = Dip Dtp'. The new fields tp and ip 
carry a U(l) charge since they transform as 

^p^^(n^p^ + n:Pn)ij^, i^.^i^^n^PL + niPn) . (2.10) 

After inserting ( p.9| ) into ( |2.7| ) we obtain an equivalent form of the action in terms of the ip 
and Ip fields 



Suy,^ = -^Yl {i'xlfiiU^xPh + U*^PR)ip^+f, - ^^+f,l^^{U*^PL + U^^PR)ij,,'j 

^ Xfl 

-vY'^xiy^PL + vfPR)iiJ. . (2.11) 



The bare mass term in ( p.5|) has now turned into a Yukawa term with Yukawa coupling 
y. Since gauge invariance is broken in the model without V field, one may add a mass 
counterterm for the gauge bosons, which takes the form of a kinetic term for the radially 
frozen scalar field in the gauge invariant extension, cf. eq. (|2.3|) above. 



Suy = ^ E [y:U,xV.+[. + V:_,f,U;,V,] . (2.12) 



Models like ( p.ll| ) have been studied in the global symmetry limit with gauge interactions 
turned off, i. e. at = with t/^^. = 1, where they reduce to pure Yukawa models. The 
phase diagram of the staggered fermion version (ISF) of the model is displayed in fig. |T]. The 
full details of this model and its phase diagram will be explained in sects. 3 and 4. Beside 
the ferromagnetic (FM) and antiferromagnetic (AM) phases at large positive and negative 
K there are two different symmetric phases in which the magnetization v = {J2x Vx) / volume 
vanishes. These two phases in fig. |I], PMW and PMS, are separated by a phase transition 
line ?/c(k), of presumably second order. The phase diagram is invariant under 
refiection which interchanges the FM and AM broken phases, cf. sect. 4. The details of the 
phase diagram depend on the charge of the scalar field in the Yukawa coupling: instead of 
the charge two field V"^ we could use for the discussion of the symmetry restoration also the 
charge one field V . A phase diagram for fermion-Higgs models where the k — > — /t symmetry 
is absent and the fermions couple to V in stead of can be found in ref. |jl8[. There the 



two symmetric phases are separated by a funnel-like region containing a ferromagnetic phase. 
Such non-universal features are, however, not important for our purpose here. 

At small y < ydi^) the Yukawa model has three scaling regions: (A) One can approach 
the FM(W)-PMW phase transition from within the FM(W) phase. Then the spectrum 
contains massive fermions and in addition there is a Higgs particle and a massless Goldstone 
boson. (B) If one approaches the FM(W)-PMW phase transition from within the PMW 
phase the spectrum contains massless fermions, but also two light scalar particles associated 
with the V field. (C) Anywhere in the PMW phase, away from the phase boundaries, the 
scalar particles decouple and even though y > 0, the spectrum contains only free massless 
fermions. These properties agree qualitatively with treating the Yukawa coupling term in 
the action ( p.llD as a perturbation. By interpreting U^^ in ( p.llQ as an external gauge field. 



we draw the important conclusion that the fermions at small y are charged. 

From the phase diagram of the Yukawa theory at small y we infer corresponding scaling 
regions (A-C) in the full theory with dynamical gauge fields, at small gauge coupling g. In 
(A) we have the Higgs phenomenon and massive fermions, in (B) we have electrodynamics 
with massless fermions and in addition massive charged scalars, while in (C) the scalars 
decouple. In the region (C) the desired gauge symmetry restoration takes place and we 
recover the target model: massless fermions interacting with the U(l) gauge fields. 



For large y {y > ydi^)) we do not recover the physics of the target modeL Previous 
investigations have shown that the effective action which describes the low energy physics 
in the PMS phase is more like eq. (|2.7| ) and not like (|2.11| ). The fermion spectrum contains 



only the ijj' fermions whose couplings to the bosonic particles vanish as a power of the lattice 



distance [0. The fermions are neutral under the U(l) gauge symmetry. They have 
masses, even in the PMS phase, which are generically of the order of the cut-off. Depending 
on details of the lattice model, the mass in the PMS phase can be tuned to zero (as in the 
Wilson- Yukawa models for a sufficiently large value of the Wilson- Yukawa coupling). In any 
case the scaling physics in the PMS phase is not that of the gauge invariant target model 
( p.4|) at y = because the fermions are neutral. 

In section 4 we shall consider yet another way of breaking the gauge invariance in the 
lattice models, namely, by changing the normalization of 75, making the replacement 

75 «:a75, (2.13) 

with ka 7^ 1- This breaks gauge invariance on the lattice, which is compact U(l), although 
in the classical continuum limit the model is still invariant under non-compact gauge trans- 
formations. Nevertheless, we shall show that in the quantum case there is a wide range of 
values of ka for which there is a PMW phase where symmetry restoration will take place. 

Summarizing, our arguments for symmetry restoration in the chiral gauge-fermion models 
with staggered fermions depend crucially on the existence of a PMW phase in the global 
symmetry limit. This phase has the massless fermions of the classical target model, which 
interact gauge invariantly with the gauge field when it is switched on as an external field. 
In the interior of the PMW phase (region (C)) there are no scaling scalars. Hence, for 
dynamical gauge fields this may be considered as a satisfactory regularization of the target 
model. 

Strictly speaking the models studied here numerically are not chiral since the charges 
Qf are chosen, for numerical reasons, such that the models are equivalent to vector models. 
However, the existence of a PMW phase depends presumably only on the relative smallness 
of the symmetry breaking and not on the chiral properties of the theory. 



3 Staggered fermion models with axial- vector couplings 

In this section we shall introduce lattice versions of the continuum model ( p.4| ) which allows 
for a chiral set of charges qf, using staggered fermions. This is most easily done in terms of 
4x4 matrix fields and \E'^ , where the indices a and k act as Dirac and flavor indices, 
respectively Consider the following ansatz for a lattice action, 

-l/ETr{^>^}. (3.1) 

X 

The projectors Pl,r on the left- (right-) hand side of \E'' project on eigenstates of 75 in the 
Dirac (flavor) space. The t/^^, . . . , may carry different representations of the gauge 
group. We restrict ourselves here to 75 as the only non-trivial flavor matrix. 

As it stands, the action ( p.l| ) has the same local gauge invariance as the Np = 4 classical 
continuum model (|2.4| ) if all components ^'^'^ and would be independent degrees of 



freedom. However, such a model would have fermion doublers and would lead to a vector- 
like gauge model for all choices of f/^^, . . . , U^^. The fermion doublers are situated at the 
boundary of the Brillouin zone in momentum space, i.e. near momenta = ±7r. After 
restricting the momenta of the matrix fields such that 

- vr/2 <p^< +7r/2 , (3.2) 

we loose the fermion doublers of the matrix field, but also gauge invariance. It is clear, 
however, that in the classical continuum limit, where the fields and their momenta go to 
zero (in lattice units), eq. (|3.1| ) with the momentum restriction ( p.2|) goes over into the 
corresponding continuum model. Hence the full gauge invariance gets restored in this limit. 

One might think that the cut-off in momentum space has to result in a non-local action. 
However, it is possible to express the action in a form that is local, using staggered fermion 
fields Xx and Xx- The connection with the and \E'^ is made by writing 0, 14 



b 



where 7^ = 71^72^73^74'' and the sum extends over the 16 corners of a unit lattice hypercube, 
bfj_ = 0,1. The substitution of ( |3.3|) into the action (|3.1| ) leads to a local action. In momentum 
space we have the relation 



^Up) = Z{p)Y,T^^,,{p)x{p + 7r6), -7r/2 <p^< +7r/2 , (3.4) 

b 

where TaK,b = Z]c6xp(i6c7r)7^^/8 is a unitary matrix and Z{p) is a non-vanishing function 
in the restricted momentum interval. This relation is derived in the appendix. It shows 
clearly that in the restricted momentum interval the Fourier components of the matrix field 
are independent. 

By choosing appropriate link fields t/^^, . . . , Ul^^, the fermions couple in different ways 
to the gauge field. For example, with 

u^' = u^' = u;:^ = u^^ , (3.5) 

the four fermions couple vector-like; with 

= u;:"" = u, , u^"^ = = 1 (3.6) 

we have a left model in which only the four left-handed fermions couple to the gauge field; 
with 

K'' = Ui"" = U, , f/«^ = U^^ = U; (3.7) 

we have a four fermion version of the axial- vector model (|2.4| ) with all = -|-1 and finally 
with 

jjLL ^ jjRR ^ ^ ^LR ^ ^RL ^ (3 g) 

we have the axial- vector model ( p. 41) with gi 2 = +1 and 52,3 = ~1- 

The key issue of gauge symmetry restoration can be investigated in model ( p.l| , pJ|) . It 
is advantageous to choose this axial vector model, and not the left model (|3.1| , |3^) , because 
the staggered version of this model has a larger lattice symmetry group which reduces the 
number of counterterms in the scaling region. Generalizing f/^ in (|3.71 ) to actions with 

arbitrary integer charge q, it is straightforward to construct chiral anomaly free models by 
combining several actions (|3.1|J3.7| ), e.g. 5(g = +1) + 4(g = —2) + l(g = +3). For numerical 



reasons we shall however study only the simple l(g = +l) + l(g = — 1) anomaly free model, 
which is equivalent to a vector theory (eight flavor QED, for y = 0), obtained by adding 
mirror fermions (cf. sect. 4). Model ( |3.8|) with an equal number of +1 and —1 axial charges 
can be easily modified such that it is exactly gauge invariant for y = (see below). This 
gauge invariant version is the earlier mentioned ISF model which is a reduction of the naive 
fermion model in eq. ( |2.5| ) and which serves here as a reference model. 

Similarly, one can construct staggered fermion models involving arbitrary coupling of the 
Dirac-fiavor components of the staggered flavors to gauge or Higgs fields in a straightforward 
manner such that the target models are recovered in the classical continuum limit. We have 
shown that this works satisfactory in the two-dimensional U(l) version of the equal charge 
axial model [0. For the Standard Model and Grand Unified Theories like SO(IO) and 
SU(5) and the axial- vector model one can write down models in which the staggered fermion 
symmetry group is preserved ||14[. This strategy of coupling the Dirac-fiavor components of 
the staggered flavors has also recently been successfully applied to a fermion-Higgs model 



After the substitution (|3.3| ) and working out the trace in (|3.1| ) one obtains the action in 
terms of the staggered fermion fields. For the equal charge axial model ( |3.1| , |3.7| ), we find 



5i 



u,x' 



E 

Xfl 



C/ix E Vfix+b{Xx+bXx 



+b+fi Xx+b+fiXx+b) 



' '"^ pLX -. „ 

E ^5x+b{V,ix +cXx+bXx+c+il VlJtx+bXx+b+[iXx+c) 
b+c=n 



l/ExLxL, (3.9) 



where 

c^3; = Re?7^^, s^^ = Imf/^^, n = (1, 1, 1, 1) . (3.10) 

The staggered sign factors 77^^. are defined as rj^x = (— 1)^i+-- +^a'-i. The factor rj^x = 
~V4^x V3x+i V2X+3+1 '7ia;+2+3+4 = —(—1)^1+^3 represents the Dirac 75. In the classical con- 
tinuum limit this action describes Np = 4 flavors of axially coupled Dirac fermions with 
Qf = +1> / 



A. 



For the axial model ( p.l| , p.8|) with gi 2 = +1 and ^3^4 = —1, we find similarly. 



S, 



u,x' 



r) E \_'^fJ'xVfJ-x{XxXx+fi Xx+fiXx) 

+ '^^^lx£xV^J.x{XxXx+^l ~ Xx+fiXx) ~ V E XxXx 



(3.11) 



The sign factor Ex = (— i)^i+---+^4 represents the product of the two 75's in Dirac and flavor 
space. From the substitution (3J) we would find hypercubical averages as in (3J) above, 
but we have chosen to replace those according to jQj^bfx-b fx- This leads to the same 
classical continuum limit and makes the model invariant under local gauge transformations 
(see below). The form ( |3.11| ) also follows by applying the familiar 'spin diagonalization' 
transformation [l^ to the naive fermion action (|2.5|) , which produces a sum of four actions 
equivalent to (p.ll|) , and keeping only one of these. 

In the equal charge axial model (|3.9|) the projection on eigenstates of the Dirac 75 leads 
to non-local couplings in the 775 term of (|3.9|) , which prevents the gauge invariance of the 
naive fermion model ( |3.1| ) to be carried over to its staggered fermion reduction. Hence we 
shall refer to this model as the non-invariant staggered fermion (NISF) model. In a similar 
fashion we would loose gauge invariance in the left model (|3.1|, p.6|) with staggered fermions. 



In the model |3.8|) however, the projection involves an additional 75 in flavor space 
and the model is invariant under the local 'e-symmetry' of the staggered fermions: the action 
( P . 1 1|) with y = is invariant under the local flavor non-singlet U(l) transformation, 



exp{iuj^ex)x^ , x'x Xx^wi^^x^x) , U^^ exp{-iuj.j;)U^x exp{iu:r+f,) . (3.12) 



Therefore we shall call the model ( p. 11 ) in the following the invariant staggered fermion (ISF) 



model. As delineated in sect. 2, gauge symmetry is restored in this model also for y > 0, 
provided y is smaller than yc- By making the transformation of variables, x'x — ^^p{'>'^x£x)Xx, 
x'x — Xx ^^vi'i^x^x) with V = expiujx we find the Yukawa form of the action in terms of the 
X and X fields 



^U,x ~ n $Z \^lixVfix{XxXx+fi Xx+fiXx) 

Xfl 

~(" '^^fj.x^xVfJ.xiXx Xx+fi — Xx+iiXx) \ — y J2[^eV,' + ^sJ^V,']xxXx, (3.13) 



with the U(l) charge two scalar field V"^. We remark that now, in contrast to the naive 
model, the fermion measure DxDx is not invariant under the transformations (|3.12|) . The 



invariance of the measure is however recovered, after adding mirror fermion fields, which is 
required also for the use of the Hybrid Monte Carlo Algorithm (HMCA) in the numerical 
simulations, cf. sect. 4. So the final ISF model has four qj = +1 and four qf = —1. 

In contrast to the ISF model the NISF model has lost its gauge invariance. The model 
can be thought of as being (almost) gauge invariant for the low momentum components of 
the gauge field while the symmetry violation increases for the high momentum components. 

Following the procedure explained in the previous section, the gauge degrees of freedom 
in the non-invariant model can be interpreted as scalar fields and the action with those scalar 
fields is obtained by replacing the sine and cosine factor in ( |3.1CI|) by 



c^a; = ReV;*?7^^\4+^ , Sf,.^ = lmV*U^^V^+f, . (3.14) 

In our numerical work we shall furthermore restrict ourselves to the zero gauge coupling 
limit Ufj_ = 1. 

We have also introduced a further modification of the models that provides an additional 
parameter to monitor a possible restoration of gauge invariance. As can be seen in the 
staggered action (^.91), the 7^ part of the target action leads only to a one-link coupling 
among the x' and x' fields, whereas the 7/^75 part involves three-and also five-link couplings. 
Therefore we expect that the 7^75 part of the action renormalizes differently from the 7^ 
part and we may have to compensate for this with a finite renormalization factor in front of 
the 7^75 term. This leads to the replacement 

Sfix AS fix (3.15) 

in the action (|3.9| ). We can also modify the ISF model in this way. For 7^ 1 the gauge 
invariance is broken and this model may be more similar to the NISF model. 

For simplicity we use the same ka for the three and five link couplings. The appearance 
of both three and five link couplings in the 7^75 term, is an awkward feature of the action 
( p.9|) , because the three and five link terms may renormalize differently. Furthermore, a 
canonical construction of a transfer matrix requires that the couplings are confined within a 
hypercube. We can arrange for this, by replacing the five link terms in (p.9|) by equivalent 



three link terms, which have the same classical continuum limit. We have sometimes used 



such a modified form of the action instead of (|3.9| ), but did not find a significant difference 
in the measured observables. 

The crucial question now is whether the breaking of gauge invariance in the NISF model 
is sufficiently weak that we can have symmetry restoration in a PMW scaling region. The 
mechanism for this would be similar to that in the ISF model, for which the gauge degrees of 
freedom Vx can be explicitly transformed into a Yukawa term (cf. ( p.l3| )), leading to charged 
massless fermions in the symmetric phase at weak Yukawa coupling. For the NISF model 
we cannot explicitly transform the V^'s into a Yukawa coupling, but the same mechanism 
may take place effectively in the low momentum modes where the symmetry breaking is 
weak. The high momentum modes might primarily lead to a 'renormalization of 75', to be 
compensated by ka- We can get an idea of the sensitivity to a mismatch in ka values by 
studying the effect of 7^ 1 in the ISF model. A major danger to the above mentioned 
scenario is that unlike the case of the massive Yang-Mills model and the ISF model, we 
cannot easily control the strength of the gauge symmetry breaking in the NISF model. 



4 Phase diagram of the ISF model 

To determine the phase structure of the ISF and NISF models we shall compute numerically 
the dependence of a number of local order parameters on the bare coupling parameters k, y 
and ka- We use the Hybrid Monte Carlo method, which requires a positive definite fermion 
determinant. The fermion matrix for the ISF and NISF models can be obtained from the 
actions ( p.ll|) and (|3.9| ) by writing them in the form 



SisF = -x'Misfx' , Snisf = -x'Mnisfx' ■ (4.1) 

The fermion determinant, Det M, is not positive definite, but this we cure by replacing it by 
Det (M^M) which is manifestly positive definite. This amounts to adding a mirror fermion 
field to the action, which couples to the matrix M^Ar)75F- With this extra field both the 
ISF and the NISF model have equal number of fermions coupling with charge qf = +1 and 
Qf = —1, which makes them vector-like and anomaly free. 

We should mention here that the actions (|3.11|) and ( p.9|) for the ISF and NISF model 
are invariant under the discrete e-symmetry 

x'x ^xX'x ' x'x X'x^x , Kr exVx , -K , (4.2) 

which implies that the n-y phase diagrams of the two models have to be invariant for k —k. 
Furthermore the actions are invariant under the symmetry 

x'x i^xx'x , x'x Xx^^x , y -y , (4.3) 

which implies that the two phase diagrams are also symmetric around the y = axis. It is 
therefore sufficient to restrict ourselves to the half plane y > 0. 

The numerical simulations were carried out mainly on a 4*^ lattice. A few simulations 
were also performed on larger lattices (6^ and 8^). We used periodic boundary conditions 
for the scalar fields. For the fermion fields periodic boundary conditions were used in the 
spatial direction and antiperiodic boundary conditions in the time direction. We typically 
accumulated a statistics of 500-3000 HMCA trajectories at each point in the phase diagrams, 
depending on the autocorrelation time for the various observables. The step size in the 
HMCA was adjusted such that the acceptance rate was between 60 and 80%. 



The local observables which we shall consider are: 
the magnetization, 

' ■ 1 



rot 



(4.4) 



where V is the lattice volume and the subscript "rot" indicates that the standard rotation 
technique has been applied to account for the drift of the magnetization vector {1/V) J2x ^ 
on a finite system 



the staggered magnetization, 



Vst 



1 



77 E^-^^ 



V 



2: 1 rot 



(4.5) 



the energy associated with the kinetic term ( [^.121 ) for the scalar field 



= (^,T.MV:V.,,)) = ^|:ln^ , (4.6) 



AVj:^ 7 8VdK 

with Z the partition function; and 

the energy associated with the mass term for the fermions 

It also proved to be useful to keep track of the number Ncg of Conjugate Gradient Algorithm 
iterations, needed for an inversion of the fermion matrix to a given accuracy. This quantity is 
large, if the fermion matrix has eigenvalues close to zero. The observables v and Vgt are order 
parameters for ferromagnetism and antiferromagnetism and are non-zero in the FM and AM 
phases, respectively. From the K-dependence of v and Vgt we determined the position of 
FM-PM and PM-AM phase transitions shown in fig. ^ A zero fermion mass would be a 
clear signal of the PMW phase, but it would require larger lattices to measure this mass 
reliably. Instead we have read off the position of the PMW-PMS phase transition from the 
^/-dependence of the energy Ey. The energies E^^ and Ey can be used to get information 
about the order of the phase transition. 

We now discuss the numerical results for the ISF model in more detail. In fig. |^a we have 
plotted the energy Ey as a function of y for k = and several values of the 75 renormalization 
factor ka (cf. (|3.15| )). The open squares were obtained at = 1, in which case the model is 
equivalent to the Yukawa model ( p.l3| ). The energy Ey has its steepest slope at the PMW- 
PMS phase transition which for = 1 is close to y = 1.4. At this point the curvature 
changes also its sign. It appears that the phase transition is of second order, which we infer 
from the absence of a jump in Ey on larger lattices and the absence of a hysteresis effect in 
thermal cycles. The energy Ey is very small, but still non-zero in the PMW phase, increases 
at the PMW-PMS phase transition and is large in the PMS phase. The approximate fall off 
(X 1/y at large y can be understood from the eigenvalue distribution of the fermion matrix 
MisF- We can write Ey = {l/V){TTMf^p) = {l/V){Zl/{X{y) +y)). The sum extends 
over all eigenvalues A of the off-diagonal part of the fermion matrix Mjsf whose eigenvalue 
distribution turns out to be almost independent of y (see below). 

The squares in fig. ^ show the y-dependence of Nqg- This quantity exhibits a maximum 
at the critical value of y indicating that the PMW-PMS phase transition is associated with 
the occurrence of small eigenvalues of the fermion matrix Mjsf- We will come back to this 
point at the end of this section. 



As we mentioned in the previous sections, it is interesting to investigate the ISF model 
with different n^. Gauge invariance is broken for 1 and we loose the equivalence of the 
models (|3.11|) and ( p.l3|) . In the figs. ||a and b we have also included the results for Ey and 
Nqg for several values of 7^ 1. Interestingly it can be seen from the shape of the curves 
that a PMW-PMS phase transition emerges for all ka> 0. The position of the PMW-PMS 
phase transition shifts to smaller values of y when the value of ka is lowered as is obvious 
from the y-dependence of the points with the steepest slope in fig. ||a and of the maximum 
in fig. 

In fig. ^ we have plotted the approximate position of the PMW-PMS phase transition, 
as determined from figs. ^ and b for n = into a ka-V phase diagram. The data points fall 
within error bars on a straight line through the origin with slope ~ 0.77. We encounter here 
another example of restoration of gauge invariance: The gauge symmetry is broken in the 
lattice action, but we expect that the PMW phase leads to the same scaling physics for all 
ka > 0. Only at ka = 0, where the 7^75 part disappears completely from the action (|3.11 



the curve for Ey in fig. ^ shows no indication for a change of curvature at small y and the 
PMS phase extends most probably down to y = 0. 

Some additional insight about the phase structure is obtained from the eigenvalue spectra 
of the fermion matrix Mjsf- The eigenvalues were computed with the Lanczos algorithm. 
In fig. H we have plotted the eigenvalues of the off-diagonal part Mjgp = Mjsf — 1/1 in the 
complex plane, for four independent scalar field configurations which were generated with 
the HMCA on a 4^ lattice at k = 0, y = 0.5 and ka = 0.3. We found empirically that the 
eigenvalue distributions obtained at other values of y, with ka and k kept fixed, have almost 
identical shape. In particular they always cut the real axis at the same distance h from the 
origin, which is ~ 0.37 for ka = 0.3. 

The distribution of Mjsf for some non-zero value of y with ka kept fixed is then approx- 
imately obtained by shifting the distribution shown in fig. ^ by an amount y along the real 
axis. This explains the observed behavior of Ncg in fig- 0b. For y < h the origin is situated 
within the hole and Nqg is small since there are no eigenvalues near the origin. It develops 
a peak when y ^ h, because eigenvalues come close to origin and it decreases again when y 
is increased beyond h. We conclude that yc ~ h, which is consistent with the value obtained 
from fig. ^. 

The width h of the eigenvalue distribution increases for increasing ka- For = 1 we 
find the eigenvalues to be arranged along the boundary of a circle and at larger ka along 
the boundary of an ellipse with the principal axis on the real axis. At = the fermion 
matrix is anti-hermitian and the eigenvalues fall exactly on the imaginary axis. This gives 
additional support for the absence of a phase transition in y at ka = 0. 

As is illustrated with this discussion, the appearance of a PMW phase at small y is 
associated with a hole in the eigenvalue distribution, which gives rise to the concave y- 
dependence of Ey for y < yc and to the peak in the number of CG iterations at y = yc- We 
shall also look for these features in the NISF model to establish if a PMW phase emerges 
there, in which the gauge symmetry breaking due to the high momentum scalar modes is 
restored. 



5 Phase diagram of the NISF model 

We have seen in the previous section that the existence of the symmetry restoring PMW 
phase of the ISF model is stable against large deviations of ha from its preferred value 1. 
This suggests that a fine tuning of ka in the NISF model may not be necessary and it should 



suffice to choose ka sufficiently large such that the PMW phase becomes clearly visible. In 
this section we discuss the phase diagram of the NISF model, mostly at fixed k = 0, but for 
various values of ka- (Preliminary results for k 7^ were reported in the first reference in 

In fig. P we have plotted the energy Ey obtained in a thermal cycle as a function of y for 

= 1. The behavior is quite different from the one observed for this quantity in the ISF 
model. The energy exhibits a jump at ?/ ~ 0.8 and a clear hysteresis effect is observed in the 
thermal cycle, indicating the presence of a first order phase transition. A similar behavior 
is found for 7^ 1 and the position of the first order phase transition in the ka-V phase 
diagram is shown in fig. |[ Note that the phase boundary does not go through the origin as 
in fig. and that there is even a phase transition at = 0. 

The properties of the NISF model in the weak coupling region are unusual and quite 
different from those in the PMW phase of the ISF model. From the behavior of the energy 

and the order parameters v and Vgt we find that the system can go into one of the following 
different states when lowering the parameter y from the PMS phase across the first order 
phase transition: 
FM : v^O,Vst = and E^ > 0. 
PM+: V = 0,Vst = and E^ > 0. 
PMq: v = 0,Vst = and E^ = 0. 
PM_: v = 0,Vst = and E^ < 0. 
AM : v = 0,Vst^O and E^ < 0. 

We shall therefore denote the weak coupling region in the following as a multi-state (MS) 
region. Also the number of CG inversions Nqg shows a qualitatively different behavior in 
the NISF model compared to the ISF model. The quantity Ncg grows when approaching 
the PMS-MS phase transition within the PMS phase, then jumps to a small value when 
crossing this transition and remains approximately constant within the MS region. 

Fig. 1^ shows the energy E^ as a function of y at k = and ka = for several thermal 
cycles started at y = 1.0 in the PMS phase. The jumps in E,^ show that the FM-PMS, 
PM+-PMS, PM_-PMS and AM-PMS are of first order. Also the PMq-PMS phase transition 
is of first order since the energy Ey exhibits a jump for all these transitions. The appearance 
of five (meta) stable states in the PMW phase indicates that the inclusion of the fermions 
produces additional (local) minima in the effective action for the scalar field Vx- The k —k 
symmetry ( [4. 21 ) relates the PM± states and the FM with the AM state. Even on very small 
lattice volumes the local minima appear to be separated by very high energy barriers in 
configuration space since we have never observed tunneling events between the various states 
of the MS region. The frequency distribution of the states which results after a repeated 
crossing of the PMS-MS transition indicate that system has a slight preference for the PMq 
phase. 

In contrast to the PMW phase of the ISF model, the MS region in the NISF model 
persists for < y ^0.7 at ka = 0. This indicates that this region is different from the PMW 
phase in the ISF model. It is not excluded, however, that a new phase, which harbors gauge 
symmetry restoration, is still present for smaller values of y, hidden inside the MS region. 
After scanning the various states of the MS region in the y and ka directions (0 < < 2, 
0.1 < y < yc{K,A)) on an 8^ lattice we did not find any indication for the emergence of 
another phase transition. 

At Ka = the 7^75 term in the fermion action is absent, and the only difference between 
the NISF model and the ISF model is the hypercubical average over Re{V*Vx+fi) in the 
one-link 7^ part of the action. This implies that the phenomenon of the MS region and the 
first order phase transition is induced by this hypercubical average. It suggests to study also 



a modified model with the replacement 



- ^ Re(\/;-.V;-5+A) - Re(\/;\4+^) (5.1) 

in the one-link term of the action ( p.9| ). In the classical continuum limit this modified model 
leads to the same target model. The modification also makes the NISF model more similar 
to the ISF model. 



6 Phase diagram of the modified NIFS model 

From the behavior of the observables f , Vst and we obtain the phase diagram in the n-y 
plane shown in fig. H, which applies to the case = 2. In contrast to the phase diagram of 
the NISF model without the replacement ( |5.1| ), we find no MS region at small y. The phase 
diagram contains only the strong coupling phases FM(S), PMS and AM(S). The FM(S)-PMS 
and PMS-AM(S) phase transitions are most likely of first order (dashed line) for y ~ 1.4 
and of second order (full line) for y ~ 1.4. The gap in E,^, v and Vgt is seen to increase when 
y is lowered. 

The K-y phase diagram at smaller values of ka looks very similar to the one shown in 
fig. H, except that the FM(S)-PMS and PMS-AM(S) phase transitions appear to change from 
first to second order at a smaller value of y. For = it appears to be second order for all 
y. This indicates that the point where the FM(S)-PMS and PMS-AM(S) phase transitions 
change their order has a similar /tA-dependence as the phase transition line in the ISF model 
shown in fig. ^. 

The apparent absence of a phase transition to a PMW phase at small y can be illustrated 
from the ^/-dependence of Ey. In fig. ^ we have displayed the energy Ey as a function of y for 
two values of ka and k = 0, on a 4^ lattice. The curve at = (open squares) is identical 
with the one shown already in fig. 0, for which the PMS phase extends down to y = 0. The 
maximum of the curve at Kyi = 2 (filled squares) is only shifted a little bit to larger y, which 
indicates that the 7^75 term needs a large prefactor ka to become important. If we try to 
fix this ka by matching the peaks in fig. ^ to those for the ISF model shown in fig. H, we 
see that the curve at = 2 of fig. ^ may be compared with the curve at ka = 0.3 in fig. ^. 
This implies that a relative factor ha ^ 7 is needed to scale the 7/^75 term in the modified 
NISF model to the same strength as in the ISF model. There is, however, no indication 
that the curvature d'^Ey/dy'^ changes sign, as for the curve at ka = 0.3 in fig. ^. Instead 
the energy appears to vanish linearly when the value of y is reduced. This indicates that a 
phase transition does presumably also not emerge when the volume is increased. Also in the 
other observables we did not find a sign for a phase transition at non-zero y. 

These results show that by the replacement ( p.l| ) we have managed to remove the MS 
region, but unfortunately we are now left with a phase diagram with only strong coupling 
phases. 

According to the discussion at the end of sect. 3, the symmetry restoration and the 
appearance of a PMW phase is expected to depend crucially on the dynamics of the low mo- 
mentum modes. The small 4^^ lattice, however, can only accommodate a few low momentum 
modes. To improve on this, we we have included at = 2 also some results on a 6^ lattice. 
All data points lie systematically slightly below the results on the 4"^ lattice, but also in this 
case there is no indication for a sign change in the curvature at small y which would reveil 
the emergence of the PMW phase. A simulation at larger values of ka was not possible with 
our resources of computer time since the number of CG iterations increases when ka is raised 



and at the same time the statistics has to be increased to account for larger autocorrelation 
times. 

Another possibility to enhance the low momentum modes of the scalar field, is to increase 
the value of n towards a second order phase transition. So far we have considered k = 
where the scalar field correlation length is presumably of order one. Therefore the high 
momentum components in Vx are very abundant and might prevent restoration of gauge 
invariance. It is therefore important to study the phase structure in a region where the 
scalar field correlation length is larger. In the pure scalar U(l) model the correlation length 
of the scalar particles in the PMW diverges when k, y k,^^^~^^^ . Also in our model we could 
hope that by choosing n sufficiently close to Kc, we can increase the scalar field correlation 
length enough that a PMW phase opens up. As delineated above the FM(S)-PMS phase 
transition at small y is however of first order. This implies that the correlation length of the 
scalar particles stays bounded and our runs at larger k show indeed no qualitative differences 
from the results for n = 0. This is illustrated by the K-dependence of Ey. Interestingly we 
find that Ey does not depend on k for a given value of y and k < k^^'^'^'^^ . In fig. ^ we have 
monitored the K-dependence of Ey for three representative values of y. This shows that a 
^/-dependence of Ey as shown in fig. ^ for ft = will also hold for other k- values in the PMW 
phase. A simulation very close to the FM-PM phase at small y and a precise determination 
of the position of the first order transition is hindered by large hysteresis effects. Therefore 
we cannot exclude the possibility that on a much larger lattice the correlation length very 
close to Kf. increases enough to allow for the appearance of the desired PMW phase. 

Additional evidence for the scenario that there is no PMW phase emerging for k Kc 'is 
obtained from the eigenvalue distributions of the fermion matrix M^isf, modified here by the 
substitution ( |5.1| ). Fig. |n| contains the spectra of M'^/^p = Mj^isp — y^ for four independent 
scalar field configurations, generated with the HMCA at the point (/«,?/, k^) = (0,0.5,2). 
The M'^/gp spectra obtained at other values of y and n, k,^^^'^^'' < k < k^^^'^^^ with 

= 2 kept fixed, look similar, except that for y ~ 0.5 the width of the distribution on 
the real axis vanishes approximately proportionally to y. Fig. [Tl| suggests no formation of a 
hole inside the distribution which, as we explained in sect. 6, would be an indication for the 
emergence of a PMW phase with restoration of the gauge symmetry at small y. 



7 Final remarks 

We have started our discussion with examples where restoration of gauge invariance takes 
place, the massive Yang-Mills theory and gauge-fermion models, provided that the symmetry 
breaking due to the bare mass term is not too large. If the mass parameters are raised 
beyond a critical value the system is driven to another phase with scaling physics differing 
substantially from that of the original gauge invariant target model without mass terms. We 
aimed for a similar dynamical gauge symmetry restoration in our chiral staggered fermion 
models, in which the gauge symmetry is broken by the high momentum modes. A strong 
motivation for this approach is that it should work without recourse to gauge fixing. 

In our approach gauge invariance is restored in the bare action by integrating over all 
gauge degrees of freedom (the 'longitudinal components'), which emerge as dynamical scalar 
fields in our lattice action, and which we want to decouple. We have applied this approach to 
a U(l) model with axial-vector couplings, and tested it in two staggered fermion realizations, 
denoted as the NISF models. 

The aim was to find in the reduction to external gauge fields, but keeping the scalars 
dynamical, a PMW phase (weak coupling symmetric phase), with massless charged fermions 



coupled gauge invariantly to the gauge field, and scalars with masses of the order of the 
cut-off. Our numerical results indicate, however, no PMW phase in either NISF model. In 
the first NISF model we could only find a PMS phase (strong coupling symmetric phase) 
connected to large values of the bare coupling parameter and a multi-state region at small 
y where the desired gauge symmetry restoration presumably does not take place. In the 
second ('modified') NISF model our results suggest that the PMS phase extends down to 
y = 0. One reason for this negative result in the modified NISF model appears to be the 
unfortunate fact that the symmetric phase at small y is separated from the broken phase by 
a first order phase transition. Therefore the scalar field correlation length remains small in 
the entire symmetric phase and the symmetry breaking due to the high momentum scalar 
modes remains too large. 

Nevertheless, we cannot completely exclude from our numerical results the possibility 
that at large values of ka, the renormalization factor of 75, a PMW phase with the desired 
symmetry restoration could still emerge in the near vicinity of the FM-PM phase transition. 
This may require much larger lattices in order to allow for sufficient dominance of the low 
momentum modes in the scaling region. It is also possible that after adding different coun- 
terterms to the action the first order phase transition between the FM and PM phase at 
small y may be changed into a second order one. One could think here of higher derivative 
couplings for the scalar field, which suppress the high momentum modes. A frustrating 
fact is, however, that in practice our simple models fail because of the rapid increase of 
computational requirements with volume and with increasing k^. 

A way out could be the following. To make gauge symmetry restoration work we need 
good control over the strength of the symmetry breaking, which was lacking in our simple 
models. A natural way to gain control would be to reduce the scalar field fluctuations on the 
scale of the fermionic lattice spacing. We could for instance use two lattices with spacings, a' 
for the fermions and a for the bosons, with a' -C a. The bosonic field values on the a' lattice 
are obtained from the a lattice by interpolation, such that the minimum wavelength bosonic 
mode is relatively smooth on the scale of a'. The ratio a' /a would then be an important 
parameter controlling the strength of gauge symmetry breaking. In the limit a' /a —* 
this would mean treating the fermions in the continuum (called 'the desperated method' in 
ref. |1|). The detailed implementation, however, of such a program is not so easy (cf. ref. |^ 
for work on interpolating lattice gauge fields). 

As an alternative one could start from a gauge fixed continuum model and regularize 
it using the lattice and staggered fermions. The model then becomes very similar to the 
Rome proposal ^ , which uses Wilson fermions. A non-perturbative test of this gauge fixing 
approach is difficult because of technical obstacles. One can still carry out simpler tests, like 
quenched simulations in gauge-fixed background configurations. Recently we have presented 
results in this direction for the case of a U(l) model with axial-vector couplings in two di- 
mensions [|T5| . 
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Appendix 

In this appendix we give the derivation of the momentum space relation (|3.4| ). The formula 
is obtained by inserting 

lb 

in the Frourier transform of eq. (|3.3|), 



The Xp ill (|A. 1|) is the Fourier transform of Xx with momentum p = q + bn. The J2q iii 



eq. ( [A.l| ) is the sum over the restricted Brillouin zone, —7i/2 < < +7r/2. After writing 
X = y + d, with ?/ on a lattice with double lattice distance, and d running over a hyper cube, 
we obtain, 

^(p) = iT7 mm exp [-i{p - q)y] exp - q)d + iqc + ?6(c + d)'K] Y^^Xq+b^ , (A.3) 
g J/ bed 

where we have used 7^ = 1 and exp^inby) = 1. Next we note that Y.y exp[i(p — q)y] = ^5p,q 
for p, q in the restricted Brillouin zone, and use 

^exp[i(c/ + c)H7aK' = Eexp(^c?67r)7i 

d d 

= 8r„,,6, (A.4) 

to obtain ( p. 41 ) with 

= n exp(zp^/2) cos(p^/2) (A.5) 



(Note that in ref. [jT4[ a factor 1/2 is missing in the expressions for T and Z{p)). By using 



the formula for irreducible representations of a group 

Y: vUgmxiaY = T^SrsS^^s^x (A.6) 

with a normalized sum, J2g = for the case of the fundamental representation of the group 
of 32 elements ±7*^, and 

5^exp[^6(c-t/)7r] = 164,,, (A.7) 

b 

it is straightforward to check the unitarity of T. 
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Figure 1: The K-y phase diagram of the ISF model. The solid lines represent the phase 
transitions between the various phases. All phase transitions appear to be of second order. 





Figure 2: a) The energy Ey as a function of y for several values of ka and k. = in the 
ISF model, b) The quantity Ncg for the same set of bare coupling parameters. Here and in 
the following we omit the error bars whenever they are smaller than the symbol size. 




Figure 4: The eigenvalues of the matrix Mjgp = Mjsf — plotted into the complex plane. 
The graph contains the eigenvalue spectra for four independent scalar field configurations 
which were generated with the HMCA at {k, y, ka) = (0, 0.5, 0.3). The lattice size is 4^. 
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Figure 5: The energy Ey as function of y for k = and Ka = 1 in the NISF model. The 
points were obtained in a thermal cycle starting at y = 1.3 with a step size of Ay — 0.02. 
Each point was obtained by taking the average over 50 HMCA trajectories. 
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Figure 6: Phase diagram of the NISF model in the ka-V plane obtained for k = 0. The 
dashed line indicates that the phase transition between the multi-state (MS) region and the 
PMS phase is of first order. 



■ FM — ^.^.^^^^^^^ 


4^ 


- PM+ — .^_-^_J 






_ PMo ^ ^ 






- PM_ — — 






_ AM ^^-^---^-^'^^'^ 





0.2 0.4 0.6 0.8 1 

y 



Figure 7: The energy E,^ as a function of y for k = and ka = in the NISF model. The 
plot contains nine thermal cycles, which were started at y = 1 inside the PMS phase with 
a step size Ay = 0.02. The plot shows that the system can tunnel into five different states 
when lowering y across the MS-PMS phase transition. Each point in a thermal cycle results 
from averaging over 50 HMCA trajectories. 
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Figure 8: The n-y phase diagram of the modified NISF model for ka = 2. The solid lines 
represent the phase transitions of second and the dashed lines phase transitions of first order. 




Figure 10: The energy Ey as a Junction of k for three values ofy and ka = 2 in the modified 
NISF model. 
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Figure 11: The eigenvalues of the matrix M^/gp = M^isf — plotted into the complex 
plane. The plot contains the eigenvalue spectra for four independent scalar field configura- 
tions which were generated with the HMCA at {K,y,KA) — (0,0.5,2). The lattice size is 



